On the Enormous Effect of Electric Field on the Crystalline Lattice of the Conductors 

with Charge-Density Waves. 
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We discuss deformation of quasi 1-dimensional conductors with charge-density wave (CDW) under 
the electric field deforming the CDW. In case of "strong" CDW-lattice interaction the effect can 
be treated in terms of the converse piezoelectric effect with enormous piezomodulus, ~ L c /\ times 
larger than that in ionic crystals (A is the CDW wavelength, L c is the CDW coherence length, 
mm scale in the sliding state). The CDW-lattice interaction is likely to be defined by the interband 
charge transfer (rearrangement of the covalent bonds) with lattice deformation, possible in a number 
of CDW compounds. The resulting effects, observed or expected, are paving a way towards new- 
principle actuators, including nanosized ones. 

PACS numbers: 71.45.Lr, 77.65.-j, 62.25.+g 



At present, one of the promising directions of studies 
in the field of the charge-density wave (CDW) physics |l| 
concerns the elastic 0] and dilation 0, 0, 0] properties of 
quasi one-dimensional conductors with CDW. Deforma- 
tion of the CDW in itself is well studied and its nature 
is in general understood. It is known, e.g., that electric 
field approaching the threshold value, E±, or exceeding it 
gives rise to non-uniform CDW deformation ("polariza- 
tion") 0,0- In this case, compression at one contact and 
strain at the other one spread over the CDW coherence 
length (more exactly, - continuity length), L c , which can 
achieve mm scale for sliding CDW @ . More poorly stud- 
ied are the effects of the CDW deformation on the crys- 
talline lattice itself Q . Probably, this lattice-superlattice 
coupling was not searched for at first , as it is not expected 
in the simple 1-dimensional model 

[11! Ed, EH Eg. 

Historically, first the CDW-lattice interaction was re- 
vealed in the studies of elastic moduli of the quasi one di- 
mensional conductors with_CDW. Decrease of the Young 



modulus, Y, up to 4%, 



|12H and of the shear 



modulus, G, up to 30% % on the CDW depin- 

ning under electric field was found. At the same time, 
such interaction is not a common feature of all the quasi 
one-dimensional compounds. E.g., K0.3M0O3, the blue 
bronze, did not show change of the Young modulus on 
the CDW depinning: \AY/Y\ < 5 • 1CT 5 [3]. 

Different interpretations of the effects have been pro- 
posed . According to [§] , if the equilibrium CDW wave- 
length, X eq , normalized by c, the lattice period along the 
chains, for some reason depends on c, a sample deforma- 
tion draws the CDW into a deformed state. In fact, in 
the absence of the CDW phase slippage (PS), the pinned 
CDW changes its period together with the crystal if the 
latter is deformed - A/c = const. If at the same time 
X eq /c const, the crystal deformation draws the CDW 
from equilibrium. So, pinned CDW contributes to the 
sample elastic modulus. In the sliding state the internal 
strains of the CDW rapidly relax, and its elastic contri- 
bution drops out 0]. 

Consider another case. Let CDW be deformed by 



an external force. Due to the interaction of the two 
"springs" - the lattice and the CDW - the crystal will 
adjust its dimensions, so that the CDW would approach 
its equilibrium, and the resulting state will be governed 
by the minimum of total elastic energy of the host lattice 
and the CDW Q. This was experimentally observed as 
metastable length states induced by application of elec- 
tric field [H or thermocycling 0]; the states, evidently, 
reveal remnant (metastable) CDW deformation. 

As the electric-field induced deformation of the CDW is 
non-uniform, study of the lattice-CDW interaction from 
the dilation of the total crystal turns out to be compli- 
cated 0, [3] ; the basic effect of the electric field on the 
lattice should be also non-uniform. So, more fruitful, to 
our understanding, would be search and study of non- 
uniform sample expansion under electric field. This idea 
is also supported by the recently found torsional strain 
induced in the quasi one-dimensional conductor TaS3 by 
electric field |a|. 

In the present Letter we will demonstrate that under 
electric field exceeding the threshold value one half of the 
sample should expand, and the other one - compress. Ex- 
perimentally, this effect could be observed, e.g., as a shift 
of the center of a sample with fixed ends. Normalizing 
this shift by the sample length and dividing by the elec- 
tric field one will obtain a value roughly ~ L c /X times 
above the typical piezomoduli of the known ionic crys- 
tals. Numerical estimate for the typical Peierls conduc- 
tor, TaS3, gives the value > 10 4 cm/V, 4-6 orders above 
the piezomoduli of the known materials. Possible mecha- 
nisms of the lattice-CDW interaction are also considered. 
The most probable explanation, as we suppose, takes into 
account the possibility of interband charge transfer un- 
der crystal deformation. The interaction is expected to 
be pronounced in compounds with covalent links, say, in 
trichalcogenides [17|, and may be negligible in the crys- 
tals with ionic bonding, like in the blue bronzes. 

Let us define the parameter <? [4, l9j characterizing the 
"strength" of the CDW-lattice interaction. Let g + 1 be 
the coefficient between the relative changes of \ eq and 
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c in the case of longitudinal deformation of the crystal 
under an external force: SX/X = (g + l)Sc/c. Clearly, 
the simple 1-dimensional model gives g = 0: in this case 
A = 27r/(j = 2/n, where n is the linear concentration of 
electrons in the metallic state. As n should change ex 1/c, 
X will change proportionally to c. E.g., for a quarter-filled 
conducting band, n = l/(2c), and A = 4c. Only g ^ 
results in the effects of the CDW-lattice interaction. The 
experimental estimates show that \g\ could be of the order 
of unity and higher 0, 0] ■ The high value of g also follows 
from the strain-induced effects in TaS3, which point out 
the transition of the CDW into the state commensurate 
with the lattice at a certain value of strain ■ We shall 
discuss the likely origin of the high g value in the end 
of the paper. Now let us calculate the profile and the 
value of a sample deformation induced by the electric 
field taking g as a given value. 

For the CDW at the depinning field, the voltage across 
the sample can be written as 



V = E t L + 2V vt 



(1) 



where L is the sample length, E t is the bulk thresh- 
old field, and V ps is the PS voltage, equal to the chem- 
ical potential shift at the contacts [lj. This relation 
means that apart from the force needed to overcome 
the bulk pinning, the excess field 2V ps /L acts upon 
the CDW chains resulting in the force per unit length 
/ = (2e/A) • (2V ps /L). This force is compensated by the 
gradient of elastic force sY c -^-t, where Y c is the CDW 
elastic modulus, <f> - its phase gain, and s is the area 
per conducting chain. If L c > L, this results in parabolic 
"sagging" of the CDW phase, (f>(x), between the contacts 
[l9j |. With cj) = at the contacts, the maximum phase 
gain in the center of the sample is: 



0(L/2) = ^eV ps Lq 2 /(sY c ), (2) 
and the displacement of the CDW is 



Sx c (L/2) = </>/q = —eV ps Lq/( s Y c ). (3) 

47T 

The longitudinal displacement (deformation) of the sam- 
ple itself, with respect to the contacts, 6x(x), follows the 
"sagging" of the CDW with the coefficient gY c /Y L Q. 
By analogy with the converse piezoelectric effect, let us 
define the "piezomodulus" of the CDW conductor, d c , as 
the displacement of the sample center normalized by the 
electric voltage applied. 

Two limiting cases worth considering. The 1st one is 
a pure short sample, where the volume pinning in the 
relation (p} can be neglected. At voltages V < 2V ps one 
can substitute V/2 instead of V ps in the relations (|2|3|) : 
normalizing Sx(L/2) by V we obtain: 



1 geLq eg L 



8tt sY l 4Y l As 



(4) 



Substituting for TaSs g = 6 0], A = 12 A, s = 20 A 2 
Q, Y L = 380 GPa 0, one obtains for L = 2 mm d c = 
5 • 10~ 6 m/V, which by 4 6 orders exceeds the know 
values for piezoelectrics. 

In the 2nd case - the more typical one - one can ne- 
glect V ps in the relation (fTJ) , and the voltage across the 
sample is dominated by the value of E t . In this case the 
"piezomodulus" appears E t L/2V ps times less: 



1 geqVp S _ eg V p . 



Ait sY L E t 2Y L E t Xs 



(5) 



This value can be, say, 1 2 order below the above esti- 
mate of d c . Still, it is very high in comparison with the 
piezomoduli of the known materials. In both cases (rela- 
tions ((31) and ([5])) the maximum expected displacement 
of the middle equals 2d c V ps , which for V ps = 3mV [l[ 
is about 300 A. Such a displacement is quite measurable, 
e.g. by means of an atomic- force microscope [2(J. 

To understand the physical sense of the enormous d c 
for the CDW compounds, we can propose an estimate for 
the piezomodulus, di, of an ion crystal - see, e.g., 2l|: 
e?" 1 ~ e/c 2 . In essence, this is an estimate of atomic 
electric fields. If, for dcfinitcncss, we take A = 4c, the 
ratio d c /di can be presented as (Yo/Yj,)(gL/4A), where 
Y = e 2 /(c 2 s). The physical sense of Y is an estimate 
of the Young modulus (or of the tensile strength) of an 



ionic crystal 



In fact, e /c is the Coulomb force 



arising on atomic displacement by ~ c (c is the scale of 
the extrapolated deformation for which the stress would 
be equal to Yj). By the way, for c = 3 A we obtain Y = 
100 GPa, in good agreement with the value for TaS3 0. 
Thus, the ratio Yq/Yl appears to be of the order of unity, 
and the relation (j4|) gives the value of d c roughly L/X 
times larger than the piezomodulus of a ionic crystals. 
The enormous d c value is associated with the macroscopic 
coherence length - the length of the expansion of the 
CDW deformations near the contacts. 

In contrast to the conventional piezoelectric effect (and 
also to the torsional strain in the CDW systems under 
electric field Q), the longitudinal non- uniform deforma- 
tion (|4I5[) has no relation to the crystal symmetry point 
group. The electric field gives rise to a gradient of the 
crystal strain 4-(Sc/c), but not to the strain 5c/ c in itself. 
Such phenomenon is known as flexoelectric effect and is 
observed in liquid crystals, as well as in solids (23[. Its 
observation, in contrast to the piezoelectric effect, does 
not require absence of the inverse symmetry of the unit 
cell. 

The calculations above imply a large value of g. To 
estimate it, one needs to consider the effect of crystal de- 
formation upon the interatomic bonding, j / implies 
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that the Fermi vector, Icf, changes in case of a longitudi- 
nal deformation of the crystal. The X eq /c variation could 
occur due to effects of the strain on the nesting - the way 
of the most integrate superposition of the Fermi surfaces, 
taking in account their curvature The nesting deter- 
mines the g-vector value. One can expect growth of cur- 
vature of the Fermi surfaces with stretching a specimen of 
a quasi one-dimensional conductor: the resulting trans- 
verse compression would increase the interchain coupling 
(just to remind, in the absence of interchain coupling, the 
Fermi surfaces would look as two parallel surfaces). We 
do not see a way to estimate theoretically the variation of 
the longitudinal component of the q- vector with the sam- 
ple strain. Neither the sign of the variation is clear. It is 
likely, that the q deviation from 2fcp is quadratic in the 
curvature of the Fermi surfaces (i.e. in the transversal 
band width). 

A more plausible explanation of the CDW-lattice inter- 
action, as we think, is related to the transfer of electrons 
(charge disproportionation) between the valency and the 
conducting-band states. A feature of the trichalco- 
genides, MX3, is the variety of the X-X bond lengths, 
In TaS3, reducing the S-S length can transform 



17] 



2S 2 ~ into S^ - with deliverance of two valency electrons 
[ltT ]. which cross over the Ta chains, to the conducting 
band ^24]. So, any crystal deformation can result in a 
change of kp 1 and, consequently, in a chage of X eq /c. 
The strain-induced charge transfer has been discussed 
also in [2] as the mechanism providing the reduction of 
the elastic moduli on the CDW depinning in NbSe3 and 
TaS 3 . 

The above explanation can give \g\ > 1. In fact, a rela- 
tively small reduction of the X-X distance, Sr xx /r xx < 1, 
could transfer 2 electrons to the conducting band. As- 
suming \Sr xx /r xx \ ~ |<5c/c| we obtain \g\ > 1. 

As we have noted above, no reduction of the Young 
modulus has been observed in the blue bronze on the 
CDW depinning: \AY/Y\ < 5 • 1CT 5 [15]. This could be 
attributed to the ion bonding in this compound, in con- 
trast to the trichalcogenides (l3j |. The difference in the 
type of Y(E) behavior of the blue bronze and of the MX3 
compounds makes us trend to the second explanation of 
the CDW-lattice interaction. 

In conclusion, the present Letter states that the crys- 
tals of quasi 1-dimensional compounds with CDW can 
inhomogeneously deform under electric fields about the 
threshold value for the onset of the CDW sliding. The 
underlying mechanism of this approximate analogue of 
the piezoelectric effect is comprised in the CDW elastic- 
ity in combination with strain-induced transfer of elec- 
trons. The possibility of such transfer is likely to be a 
feature of the Peierls quasi one-dimensional conductors 
comprising S or Se [l7| . The expected huge non- uniform 
longitudinal deformation, as well as the recently reported 
torsional strain open prospects for developing of new- 
principle electric actuators covering the areas of micro- 



and nanometer dimensions. 
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